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基于 尺度 不 变 特征 点 匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 


摘要 : 不 同 太阳 观测 站 间 存 在 观测 
Fer 有 必要 找到 不 
匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 。 


和 定位 。 开 展 了 不 同时 段 磁场 图 像 的 配 准 和 
差 等 参数 进行 了 定量 分 析 ， 实 验 结果 表明 ， 该 方法 可 以 自动 、 准 确 
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0012) 


MER! 


上 对 间 、 观 测 范 围 、 观 测 设备 等 差异 ， 为 了 方便 科学 研 


两 幅 图 像 特 征 点 进行 粗 配 准 和 粗 定位 ; 最 后 


定位 实验 ， 对 


像 的 配 准 与 定位 。 
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太阳 磁场 在 太阳 活动 中 扮演 了 极为 关键 的 


发 活动 ， 须 监测 从 太阳 光 球 到 
有 多 个 太阳 磁场 观测 设备 。 如 国 
的 口径 35cm 太阳 磁场 望远镜 (Solar Magnetic 
天 文 台 (Solar Dynamics Observatory, SDO) 的 
Imager, HMI) B 等 。 不 同 太阳 观测 站 间 存 在 观测 时 间 、 观 测 范围 
提高 太阳 物理 数据 使 用 效率 ,有 必要 找到 不 同 观测 
法 。 另 外 ， 当 两 观测 设备 的 观测 范围 不 同时 ， 还 需要 实现 “小 图 
的 太阳 磁场 望远镜 得 到 的 是 局 部 太阳 磁场 图 


便 科 学 研究 ， 


定位 。 例如， 


怀柔 太阳 观测 基地 上 


动力 学 天 文 台 的 日 震 及 磁场 成 像 
后 者 ， 这 就 需要 找到 一 种 能 够 将 
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晃 的 磁场 和 流 场 结构 及 
内 怀柔 太阳 观测 


仪 得 到 的 是 全 


全 日 面 图 像 


匹配 点 对 


E 提 取 


数 、 匹 配 ; 


同 观测 站 间 太 阳 图 像 自动 配 准 与 定位 的 方法 。 提出 了 基于 尺度 不 变 特征 点 
先 , 对 原始 图 像 进行 对 比 度 增 强 、 
采用 尺度 不 变 特征 检测 算法 提取 两 图 的 尺度 不 变 特征 点 ; 然后 ， 寻 找 同名 点 的 方法 对 


降 采 样 等 预 处 理 ; 


EMK, 


， 基 于 粗 定位 区 域 实 现 太阳 磁场 图 像 的 精确 配 准 


、 人 快速 地 实现 太阳 磁场 图 


已 经 建 


匹配 误 


角色 。 为 了 预测 耀 班 和 日 晃 物 质 抛射 这 样 的 爆 
此 演化 过 程 中 。 目 前 ， 国 内 外 


站 间 太 阳 图 
像 ” 在 “大 


图 像 ”中 的 
像 ， 然 而 ， 太 阳 


基地 (Huairou Solar Observing Station,HSOS) 
Field Telescope, SMFT) 路、 国外 太阳 动力 学 
震 及 磁场 成 像 仪 (Helioseismic and Magnetic 
、 观 测 设备 等 差异 。 为 了 方 
像 自动 配 准 


的 方 


动 


日 面 太阳 磁场 图 


和 象 。 前 者 的 观测 范 目 
快速 定位 的 方法 。 


局 部 太阳 图 像 


学 院 太阳 活动 


点 实验 室 


家 自然 科学 天 文联 合 基金 (U1531247，U1331113); 科技 部 基础 工作 专项 重点 项 目 


放 课 题 (KLSA201617) 资 助 . 


目 远 远 小 于 
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太阳 


对 同一 观测 站 的 序列 图 像 进行 的 ， 其 成 像 设备 位 置 


目前 


已 有 一 些 关于 太阳 图 像 配 准 方面 的 研究 工作 ， 如 文 [外 提出 了 
图 像 配 准 算法 ， 文 [5] 采 用 基于 信息 粹 与 SIFT 算法 对 太阳 图 像 进行 了 配 准 。 这 些 都 是 针 


基于 修正 矩 的 亚 像 素 


定 ， 图 像 内 容 在 短 时 间 内 变化 较 小 ,位 


El 
移 、 旋 转 和 仿 射 变换 较 少 。 不同 太阳 观测 站 的 太阳 图 像 由 于 在 观测 时 间 、 观 测 范围 、 观 测 设 


等 方面 存在 较 大 差异 , 需要 适应 力 更 强 的 图 像 配 准 方法 进行 配 准 。 已 有 的 图 像 配 准 方法 很 
于 灰 度 的 配 准 方 法 IT、 基 于 特征 的 配 准 方法 B31 等 等 。 基 于 灰 度 的 配 准 方法 具有 精 


尺度 不 变 特征 点 匹配 的 图 像 配 准 方法 ao20 计算 量 小 ， 


较 好 
变 特 


像 在 太阳 动力 学 天 文 台 日 


征 点 


匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 , 开展 了 怀柔 太阳 观测 基地 的 局 部 太阳 磁场 图 


并 做 了 检测 实验 证 明 方 法 的 有 效 性 。 


1. 算法 原理 


度 不 变 特征 提取 算法 


4 个 步 又。 


尺度 不 变 特征 点 匹配 算法 对 局 部 太阳 图 像 在 全 日 
的 特征 点 检测 、 特 征 点 初步 匹配 、 特 征 点 精确 匹配 、 图 像 变换 矩阵 求解 


点 ， 但 计算 复杂 度 高 ， 且 对 目标 的 旋转 、 形 变 、 遮 挡 以 及 灰 度 变化 等 很 敏感 。 基 于 
对 灰 度 的 变化 、 图 像 形变 及 遮掩 等 有 
的 适应 能 力 。 太 阳 磁 场 图像 含 有 许多 内 部 结构 ， 特 征 点 较 明显 ,本 文 提出 了 基于 尺度 不 


震 及 磁场 成 像 仪 的 全 日 面 太 阳 磁 场 图 像 中 自动 配 准 与 定位 的 研究 ， 


太阳 中 自动 配 准 与 定位 , 主要 包括 尺 


1.1 SIFT 特征 点 检测 
1.1.1 尺度 空间 极 值 检测 


变换 


R 度 
Be 


一 幅 


空间 理论 是 检测 不 变 特征 的 基础 。 文 [12] 证 明了 高 斯 卷 积 核 是 实现 尺度 变换 的 唯一 


二 维 图 像 在 不 同 尺度 下 的 尺度 空间 可 以 表示 为 图 像 与 高 斯 核 卷 积 : 


L(xX,y,0)= G(x, y,o)* I(x, y) $ (1) 


其 中 ，Ccoy) 为 图 像 点 的 像素 坐标 ，7Ccy) 为 图 像 灰 度 值 ，CG(x, y, o) 为 高 斯 核 函数 。 


像 的 


E 
x+y 


1 y 
G(x, y,o) = err, 
2 2x0” 


尺度 空间 。 


为 


为 高 效 地 在 尺度 空间 检测 稳定 的 极 值 ， 文 [10] 使 


断 依 


(2) 


其 中 , o 为 尺度 空间 因子 , 它 是 高 斯 正 态 分 布 的 方差 , 反映 了 图 像 被 平滑 的 程度 ; Lx, y,o) 


尺度 空间 中 差分 高 斯 DoG 的 极 值 作 


Ho DoG 算 子 定义 如 下 : 
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D(x, y, o) =[G(x, y,ko) - G(x, y,o)]* I(x, y) 
= L(x, y,ko) — L(x, y,o) 


其 中 , 上 为 相 邻 两 个 尺度 空间 的 比例 因 


子 。 


DoG 局 部 极 值 的 检测 是 将 每 个 像素 的 DoG 值 跟 同 一 尺度 的 周围 邻 域 8 个 像素 和 相 邻 尺 


度 对 应 位 置 的 周围 邻 域 9x2 个 像素 总 


G) 


k 26 个 像素 进行 比较 。 仅 当 被 检测 点 的 DoG EKF 


26 个 像素 点 或 小 于 26 个 像素 点 时 ， 才 将 其 判断 为 极 值 点 并 保存 ， 以 进行 后 续 计 算 。 


11.2 确定 关键 点 位 置 及 尺度 


通过 拟 合 三 维 二 次 函数 以 精确 确定 关键 点 的 位 置 和 尺度 , 从 而 可 


[得 到 原 图 像 尺 度 不 变 特 


征 提 取 算 法 的 候选 点 集合 。 其 中 ， 对比度 低 的 点 对 噪声 比较 敏感 ， 而 位 于 边缘 上 的 点 难以 准 


确定 位 。 为 了 确保 尺度 不 变 特征 提取 算法 特征 点 的 稳定 性 , 须 从 


点 集合 将 这 两 类 点 剔除 。 


1.1.3 关键 点 方向 确定 


旋转 不 变性 。 关 键 点 的 方向 利用 邻 域 像素 的 梯度 分 布 特 怕 


L(x,y) 的 梯度 方向 0(x,y) 可 以 由 下 式 计算 00， 


确定 。 对 4 


尺度 不 变 特征 提取 算法 候选 


通过 确定 关键 点 的 方向 , 可 以 使 特征 描述 符 以 与 方向 相关 的 方式 构造 ， 从 而 使 算 子 具有 


O(x, y) = arctan 


1.1.4 特征 向 量 生成 


为 了 确保 描述 符 的 旋转 不 变性 ， 


后 将 平 


可 能 分 布 的 360° 梯度 方向 划 为 8 个 方向 范围 ，4 


每 个 方向 范 


围 各 包含 45" 。 在 子 区 


每 个 高 斯 图 像 ， 每 个 点 


| L(x, y+ - L(x, y-D m 
| L(x+1, y)- L(x-1, y) | 
先 需 旋转 坐标 ,保持 坐标 方向 与 极 值 点 方向 一 致 ;， 然 


域 为 4x4 的 情况 下 ， 共 计 有 4x4x8=128 个 数据 ， 即 生成 128 维 描述 符 。 方 向 与 信息 联合 不 
旋转 和 变形 的 影响 510。 


仅 增 强 了 抗 噪 性 ， 还 消除 了 尺度 变化 、 
1.2 特征 点 初步 匹配 


完成 特征 点 检测 后 ， 可 利 


寻找 同名 点 的 方法 对 两 幅 图 


像 特征 点 进行 粗 


式 距离 d 作为 图 像 间 的 相似 性 度量 以 确定 同名 点 。 欧 式 距离 二 维 表示 为 


d= Ja 


两 个 值 中 小 的 值 作为 分 子 ， 大 的 值 作为 分 母 
应 点 的 同名 点 关系 ， 称 为 匹配 点 。 阔 值 的 选 
匹配 点 越 多 ， 可 能 的 误 匹 配 也 越 多 ， 阔 值 过 小 时 ， 


0.4~0.6 之 间 最 佳品。 


i) | (y, yo) : 


1.3 特征 点 精确 匹配 及 变换 矩阵 求解 


PERES, WARK 


(5) 


先 在 待 配 准 图 像 中 找 出 两 个 与 参考 图 像 中 对 应 某 极 值 点 的 欧 氏 距离 最 短 的 点 ; 然后 以 
， 若 该 分 数值 在 给 定 的 闵 值 范围 内 ,那么 建立 对 
取 将 直接 影响 同名 点 的 数量 。 阔 
匹配 点 会 很 少 。 


值 越 大 ， 获 得 的 
大 量 实验 表明 ， 闭 值 在 
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在 完成 特征 点 粗 匹 配 后 ， 可 利用 随机 抽样 一 致 性 算法 外 消除 错误 的 匹配 点 ， 从 而 提高 配 
准 精 度 。 随 机 抽样 一 致 性 算法 ， 可 以 在 一 组 包含 “外 点 ”的 数据 集中 ， 采 用 不 断 迭 代 的 方法 ， 
寻找 最 优 参 数 模型 ， 不 符合 最 优 模型 的 点 ， 被 定义 为 “外 点 "。 该 算法 最 早 由 文 [9] 提 出 四， 在 
图 像 配 准 以 及 拼接 上 得 到 广泛 应 用 。 


一 般 地 ， 两 幅 图 像 间 的 变换 可 以 由 变换 矩阵 表示 031: 
x h h h fx 
y|=| hh hs yy ， 
1 h h, 141 
(6) 


Hp, Cx, yO 和 (x, y) 分 别 为 两 幅 图 像 中 对 应 点 的 像素 坐标 ; hi Ci=0,1,2...7) 为 变换 


和 矩阵 的 8 个 未 知 系数 。 

随机 抽样 一 致 性 算法 可 对 特征 点 进行 精确 匹配 并 求 变 换 矩 阵 。 首 先 随机 选择 4 组 对 应 点 
组 成 一 个 随机 样本 并 计算 变换 矩阵 五; 对 假设 的 每 组 对 应 计算 距离 4; 计算 与 有 一 致 的 内 点 
数 ; 选择 有 最 多 内 点 数 的 态 ， 在 数目 相等 时 ， 选 择 内 点 标准 方差 最 小 的 那个 解 (3。 最 后 根 
据 图 像 间 变换 矩阵 H, 可 以 计算 出 局 部 太阳 图 像 的 4 个 顶点 及 中 心 点 在 全 日 面 图 像 中 的 像素 
坐标 ， 从 而 完成 局 部 太阳 图 像 在 全 日 面 太阳 图 像 中 的 配 准 与 定位 。 


> 


2. 算法 设计 及 实验 分 析 


2.1 算法 设计 
基于 尺度 不 变 特征 匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 流程 如 下 : 

(1) 对 图 像 进行 预 处 理 ， 包 括 降 采样 、 黑 白 反 相 、 对 比 度 增 强 等 ， 对 图 像 降 采样 到 
1024*1024 左右 ; 
D 通过 尺度 不 变 特征 提取 算法 提取 图 像 特征 点 ; 
(3) 采用 寻找 同名 点 的 方法 对 两 幅 图 像 特征 点 进行 粗 匹配 ; 
(4) 利用 随机 抽样 一 致 性 算法 去 除 错误 的 匹配 点 ， 并 估计 两 幅 图 像 之 间 的 最 佳 变 换算 


阵 H; 


(5) 根据 仿 射 变换 和 矩阵， 实现 局 部 太阳 图 像 在 全 日 面 太 阳 图 像 中 的 初步 定位 ; 
(6) 利用 初步 定位 的 图 像 区域 和 怀柔 局 部 图 像 进行 精确 配 准 。 
算法 的 具体 流程 如 图 1。 


= 


图 像 特征 特征 变换 初步 精确 
oe 点 检 点 配 矩阵 定位 EO 定位 
测 准 


= 
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1 太阳 图 


像 配 准 与 定位 流程 


图 


Fig.l Flow chart of solar image registration and positioning 


2.2 实验 结果 与 分 析 
2.2.1 不 同 观测 站 太阳 磁场 图 像 的 配 准 与 定位 
以 怀柔 太阳 观测 基地 的 局 部 太阳 磁场 


匹配 与 定位 


给 出 


展 


开 实验 。 图 2 给 出 了 原始 的 全 日 


图 像 在 SDO/HMI 的 全 日 面 太 阳 磁 场 图 像 中 的 自动 
面 太阳 图 像 及 局 部 太阳 图 像 。 其 中 ， 
日 震 及 磁场 成 像 仪 中 得 到 的 全 日 面 太阳 磁场 图 像 , 拍摄 时 间 2012 年 11 月 12 日 05 时 48 


图 2(a) 


分 00 秒 。 图 2(b) 给 出 35cm 太阳 磁场 望远镜 外 得 到 的 局 部 太阳 磁场 图 像 ， 拍 摄 时 间 2012 


围 和 分 辨 率 有 很 大 差异 。 怀 柔 太 阳 观 测 基地 的 局 部 太阳 图 像 仅 仅 是 全 日 


的 部 分 。 全 图 像 的 尺度 830x992。 


年 11 月 12 


面 图 像 的 尺度 是 4096x4096， 怀 柔 站 局 部 


日 05 时 55 分 34 秒 。 两 者 的 拍摄 时 间 非 常 


榜 近 。 值 得 指出 的 是 ， 


两 者 的 观测 范 
太阳 图 像 中 很 小 一 


(a) 全 日 


太阳 磁场 图 像 ; 


一 
T 


2 原始 的 太阳 磁场 图 像 


) 


局 部 太阳 磁场 


=> 


像 


Fig. 2 Original solar magnetic field image(Full-disk solar magnetic field image in(a)and local solar magnetic 


两 幅 图 
准 速度 ， 先 做 降 采 村 
成 像 仪 的 全 
过 
增强 后 的 局 


图 


图 2 中 两 幅 


面 太阳 磁场 图 


降 采 样 、 


p 


部 
图 


反 相 及 对 比 度 增强 处 理 后 


太阳 图 像 如 图 


field 


图 像 的 对 比 度 不 够 。 在 进行 后 续 的 图 
像 进行 预 处 理 。 另 外 ， 由 于 


处理， 将 尺度 从 4096x4096 降低 到 
像 的 正 负 极 性 与 怀柔 观测 站 的 相 
的 全 日 面 太阳 磁场 


image (in(b))) 


震 及 磁场 成 像 仪 的 全 日 


像 配 准 前 ， 有 必要 通过 对 比 度 增强 等 对 


j 图 像 分 辨 率 太 高 ， 为 了 提高 配 


1024x1024。 最 后 ， 部 分 


震 及 磁场 
反 。 因 此 要 做 黑白 反 相 处 理 。 经 
图 像 如 图 3 (a) 。 经 过 对 比 度 


3 (b) 。 对 比 图 2 (b) 和 图 3 (b) 可 以 发 现 ， 预 处 理 后 的 太阳 


BERRE 


像 更 加 清晰 ， 


这 有 利于 后 续 的 特征 提取 与 图 像 配 准 。 
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+ 


(a) (b) 


(a) 全 太阳 磁场 图 像 ; (0) 局 部 太阳 磁场 图 像 
图 3 经 过 预 处 理 后 的 太阳 磁场 图 像 


Fig.3 Full-disk solar magnetic field image(in(a))and local solar magnetic field image (in(b))after pre-treatment 


随后 , 采用 尺度 不 变 特征 提取 算法 提取 图 像 特 征 点 。 图 4 给 出 了 经 过 预 处 理 后 的 局 部 太 
阳 图 像 的 尺度 不 变 特征 检测 结果 。 一 共 检 测 到 了 204 个 尺度 不 变 特征 点 。 由 图 可 见 ， 这些 特 


4 


征 点 较 好 地 代表 了 局 部 太阳 图 像 的 特征 。 


200 400 600 800 


4 经 过 预 处 理 后 的 局 部 太阳 图 像 的 尺度 不 变 特征 检测 


Fig.4 Scale-invariant feature detection of localized Solar image after pre-treatment 


在 对 图 3 中 的 两 幅 图 像 进行 尺度 不 变 特征 检测 后 ， 需 要 将 两 幅 图 的 特征 点 配 准 。 图 3 
中 的 全 日 面 图 像 一 共 检 测 到 了 3447 个 尺度 不 变 特征 点 ， 局 部 太阳 图 像 一 共 检 测 到 了 204 个 
尺度 不 变 特征 点 。 接 下 来 , 采用 寻找 同名 点 的 方法 以 及 随机 抽样 一 致 性 算法 找 出 两 幅 图 中 相 
互 配 准 的 特征 点 对 。 如 图 5, 一 共 找到 了 20 个 相互 配 准 的 特征 点 对 (图 中 用 红色 直线 连接 ) 。 
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根据 这 些 配 准 特征 点 对 ， 可 
-0.13 

H, =| -0.05 

0 
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初步 估计 出 两 幅 图 像 之 间 的 变换 矩阵 D, BARUTA: 


75 0.0583 727.0 
73 -0.1303 675.9) . (7) 
0 1 


根据 变换 矩阵 已 可 以 计算 出 局 部 图 像 的 4 个 顶点 及 中 心 点 在 全 日 面 图 像 中 的 像素 位 置 ， 
从 而 完成 了 局 部 图 像 在 全 日 


面 图 像 中 的 粗 定位 〈 图 5) 。 


图 5 局 部 图 像 在 全 日 面 图 像 中 的 配 准 与 定位 


Fig.5 Registration and positioning of local image in full-disk solar image 


对 日 震 及 磁场 成 像 仪 全 


日 面 图 像 进 行 降 采 样 处 理 可 以 大 大 提高 计算 速度 , 但 同时 也 会 降 


像 中 的 精确 范围 。 精 确 匹 配 
3 红色 直线 连接 ) 。 根 据 这 
体 如 (8) 式 : 


-0.5425 0.2270 547.2349 | 


氏 匹 配 与 定位 精度 。 为 此 , 将 日 震 及 磁场 成 像 仪 全 日 面 原始 图 像 中 绿色 四 边 形 对 应 的 最 小 外 
BRET Xk COLES) 与 怀柔 局 部 太阳 图 像 做 精确 匹配 ， 即 可 确定 局 部 太阳 图 像 在 全 日 面 图 


后 的 结果 如 图 6， 一 共 找 到 了 25 个 相互 配 准 的 特征 点 对 《图 中 
些 配 准 特征 点 对 ， 可 估计 两 幅 图 像 之 间 的 最 佳 变换 矩阵 H2 F 


H, =| -0.2221 -0.5300 662.1596] . (8) 


2 


0 


0 1 


由 变换 矩阵 AD 可 以 看 出 ， 两 图 像 间 主 要 存在 缩放 、 旋 转 、 平 移 等 几何 变换 ， 不 存在 透 


射 变换 , 且 可 以 计算 局 部 图 


RAY 4 个 顶点 在 全 日 面 图 像 中 的 精确 位 置 分 别 为 (2909,2704) ~ 


(2371,2484) 、(2560,2045) ~ (3097,2265) (图 6 中 用 绿色 四 边 形 标 出 ) ， 从 而 完成 了 


局 部 图 像 在 全 日 面 图 像 中 的 


精确 定位 。 另外, 还 可 以 计算 局 部 太阳 图 像 相对 于 全 日 面 图 像 的 


旋转 角度 约 为 157.9”“【〔 顺 时 针 〉; 水 平方 向 放大 比例 为 1.71 、 垂 直方 向 放大 比例 为 1.74。 
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图 6 


像 在 全 


Fig.6 Accurate registration and positioning of local image in full-disk solar image 


nkr 


BRRR (A 7 (b) ) 。 


图 7 给 出 了 局 部 放大 后 的 全 日 


图 可 见 ， 两 图 的 形状 、 角 度 等 基本 一 致 ， 但 细节 方 


和 


。 两 图 形状 、 角 度 等 基本 


就 


是 由 不 同 观测 站 得 到 的 ， 在 观 涡 


时 间 、 望 远 镜 参 数 等 方 


FE 明了 本 文 算法 的 准确 性 。 两 图 细节 方面 略 有 差异 是 


中 的 精确 配 准 与 定位 


MAPA CRT (Ca) ) 与 经 过 变换 矩阵 H 变换 
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后 的 局 


异 不 会 影响 局 部 太阳 图 像 在 全 


( 


图 像 中 的 配 准 与 定位 。 


a) 局 部 放大 后 的 全 


太阳 磁场 图 像 ; 


47 全 日 面 太阳 磁场 


= 


四 像 与 经 过 变换 后 的 


为 两 
j 存 在 差异 所 致 。 但 这 些 细微 的 差 


Fig.7  Full-disk solar magnetic field image ( in( a ) )and transformed local solar magnetic field image( in ( b ) ) 


2.2.2 不 同时 间 段 的 匹配 结果 


为 了 验证 算法 的 有 效 性 ， 对 71 组 不 同时 段 的 太阳 
对 数 与 正确 配 准 点 对 数 ， 结 果 如 图 8。 其中， 红色 实 线 圆 
点 为 正确 配 准 的 点 对 数 ， 蓝 色 虚 线 五 角 星 为 配 准 的 正确 率 ( 即 1 


E， 并 给 出 相应 的 配 准点 
， 绿 色 虚 线 


图 像 进 行 配 ; 


点 


准 的 点 对 数 与 实际 配 准 
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的 点 对 数 之 比 )， 全 日 
配 准 的 点 对 数 与 实际 的 配 准点 对 数 相等 ， 


太阳 图 


BA 
大 气 消 流 等 干扰 下 ， 怀 柔 局 部 大 


=} 


震 及 磁场 成 像 仪 全 


73 


而 导致 配 准点 对 数 和 配 i 


怀柔 局 部 太阳 图 像 较为 模糊 时 ， 配 准点 对 数 与 配 准 准确 率 较 低 。 


面 图 像 分 辩 率 均 为 1024x1024。 由 图 8 可 以 看 出 ， 大 部 分 情况 下 正确 
即 配 准 准确 率 为 100%; 少数 情况 下 正确 配 准 的 点 
对 数 小 于 实际 的 配 准 点 对 数 ， 即 存在 一 定 的 错误 配 准 。 出 现 错误 配 准 的 原因 在 于 : 怀柔 局 
面 太 阳 图 像 ， 一 个 是 地 面 拍摄 的 ， 一 个 是 空间 拍摄 的 。 
日 图 像 会 变 得 模糊 ，SIFT 算法 能 检测 到 的 特征 点 减少 ， 
准确 率 的 下 降 。 图 9 给 出 了 局 部 太阳 图 像 较 模 糊 时 的 配 准 情况 ， 
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图 8 不 同时 间 太 阳 磁 场 图 像 的 配 准 情 况 


Fig.8 The registration of the solar magnetic field image at different time 


200 400 600 800 1000 1200 1400 1600 


图 9 局 部 太阳 磁场 图 像 较 模糊 时 的 配 准 情况 
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Fig.9 Registration results of a blurry local solar magnetic field image 


2.2.3 仿 射 变换 对 太阳 磁场 图 像 配 准 结果 的 影响 


不 同 观测 站 的 太阳 磁场 图 
较 明 显 的 仿 射 变换 。 本 节 研 究 仿 射 变换 对 太阳 磁场 图 像 配 ; 


像 间 除了 常见 的 缩放 、 旋 转 和 3 


F 移 变换 , 很 多 情况 下 还 存在 比 


住 结果 的 影响 ， 通 过 定量 分 析 并 计 


算 配 准 实验 的 误差 。 从 原始 日 震 及 磁场 成 像 仪 全 日 面 图 像 中 (图 2(a) ) 抠 出 一 块 矩 形 子 图 像 。 


然后 ， 将 矩形 子 图 像 进行 


n== 


7 BE FU A E o 


仿 射 变换 包含 平移 、 旋 转 、 缩 放 、 错 切 等 变换 。 图 
配 准 误差 (单位 : 像素 ) 随 错 切 系数 和 旋转 角度 的 变化 情况 。 由 
错 切 系数 增 大 到 一 定 程度 如 超过 0.3 后 , 配 准 误差 将 迅速 增 大 ; 
较 小 。 图 11 给 出 了 旋转 角度 为 45 度 的 情况 下 , 配 准 误差 随 错 切 系 
图 11 可 以 看 出 ， 缩 小 或 放大 2 倍 以 内 ， 缩 放 对 配 准 误差 的 影 


较 小 时 , 配 准 误差 较 小 , 而 当 


旋转 角度 对 配 准 结果 影响 


数 和 缩放 系数 的 变化 情况 。 


定 的 仿 射 变换 ， 最 后 ， 将 仿 射 变换 后 的 子 图 像 来 与 降 采 样 后 的 全 


图 像 (图 3(a)) 做 配 准 ， 将 变换 矩阵 计算 得 到 的 子 图 像 在 


a 


全 日 面 图 像 中 的 位 置 与 子 图 像 


在 全 日 面 图 像 中 的 已 知 位 置 相 减 ， 即 可 准确 地 计算 配 准 ; 


和 对。 配 准 误差 7 定义 为 


DG = Xo)? +(y; = YoY 


4 


HH, (xo yio) 和 (%, y ) 分 别 为 子 图 像 的 4 个 顶点 在 全 日 


(9) 


图 像 中 的 已 知 位 置 和 经 过 变换 


10 给 出 了 缩放 为 0.6 倍 的 情况 下 ， 


图 10 可 以 看 出 ， 当 错 切 系 数 


特征 点 的 检测 与 配 准 千 
变换 等 均 具 有 一 定 的 鲁 棒 


响 较 小 ,超出 此 范围 后 配 ; 


E 误 差 将 明显 增 大 。 这 是 因为 缩小 变换 需要 丢弃 一 些 图 像 信息 ,而 
放大 变换 则 需要 插值 来 增 大 图 像 。 当 缩小 导致 丢弃 的 信息 太 多 或 放大 插值 的 信息 太 多 时 会 对 
影响 。 综 上 所 述 ， 本 文采 用 的 SIFT 配 准 算法 对 旋转 、 缩 放 、 错 切 
性 。 为 了 更 好 地 反映 仿 射 变换 对 配 准 精度 的 影响 ,本 节 的 结果 都 是 


通过 初步 定位 得 到 的 〈 对 应 算法 步骤 ave) ， 通 过 精确 定位 可 以 大 大 减 小 配 准 误差 。 
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图 10 配 准 误差 随 错 切 系数 和 旋转 角度 的 变化 


Fig.10 The change of registration error varies with the shear factor and the rotation angle 
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图 11 配 准 误差 随 错 切 系数 和 缩放 系数 的 变化 


Fig.11 The change of registration error varies with the shear factor and zoom factor 


T 


T= 


T 


采用 基于 尺度 不 变 特征 点 匹配 的 方法 对 怀柔 观测 站 的 局 部 太阳 磁场 图 像 在 太阳 动力 学 
人 台 日 震 及 磁场 成 像 仪 的 全 日 面 图 像 中 的 配 准 与 定位 问题 进行 了 研究 , 结果 表明 , 该 算法 
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对 光照 、 旋 转 、 尺 度 缩放 等 有 很 强 的 鲁 棒 性 ， 对 不 同 观测 站 得 到 的 具有 不 同 观测 条 件 、 不 同 
观测 尺度 等 的 太阳 磁场 图 像 的 配 准 有 效 且 提高 了 太阳 物理 数据 使 用 效率 。 
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Registration and Location Method of Solar Magnetic Field 


Images Based on Scale-invariant Feature Point Matching 
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Abstract: There are differences among different solar observation stations in terms of observation 
time, scope and equipment. For the convenience of scientific research, it is necessary to work out a 
way to automatically match and locate the solar images of different observation stations. A method 
to match and locate solar magnetic field image based on scale-invariant feature point matching is 
put forward. Firstly, pre-treatment for the original image is carried out by contrast enhancing and 
down sampling; Secondly, the feature scale-invariant feature points are extracted using 
scale-invariant feature point detection algorithm; Next, the feature points of the two images are 
matched and located roughly by seeking their homologous points; Finally, ,Based on coarse 
located area ,the solar magnetic image can be matchend and located exactly . The registration and 
positioning experiments of the solar magnetic field image in different time periods is carried out, 
and number of matching point pairs,matching logarithm, matching accuracy and matching error 
are defined to quantitatively analyze the experimental results. As proved by experiment results, the 
presented method can match and locate the solar magnetic images automatically, accurately and 
quickly. 

Keywords:Solar image; Image registration; Automatic positioning; Scale-invariant; Feature 


extraction 


